Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 September 2009 (MN LSJgC style file v2.2) 



Simulating the effect of AGN feedback on the metal enrichment of 
galaxy clusters 

O D- Fabjan 12 ' 3 , S. Borgani 1 ' 2 ' 3 , L. Tornatore 1 ' 2 ' 3 , A. Saro 1 ' 2 3 , G. Murante 4 & K. Dolag 5 

, 1 Dipartimento di Astronomia dell 'University di Trieste, via Tiepolo 11, 1-34131 Trieste, Italy (fabjan,borgani,tornatore, saro@oats.inaf.it) 
• 2 INAF — Istituto Nazionale di Astrofisica, via Tiepolo 11, 1-34131 Trieste, Italy 
3 INFN - Istituto Nazionale di Fisica Nucleare, Trieste, Italy 
~ , 4 INAF — Istituto Nazionale di Astrofisica - Osservatorio Astronomico di Torino, Str. Osservatorio 25, 1-10025, Pino Torinese, Torino, Italy (murante@oato.inaf.it) 
^ -* 5 Max-Planck-lnstitut fur Astrophysik, Karl-Schwarzschild Strasse 1, Garching bei Miinchen, Germany (kdolag@mpa-garching.mpg.de) 

m ■ 



O 

U 



3 September 2009 



Oh 



ABSTRACT 

We present a study of the effect of AGN feedback on metal enrichment and thermal proper- 
ties of the intracluster medium (ICM) in hydrodynamical simulations of galaxy clusters. The 
simulations are performed using a version of the TreePM-SPH GADGET- 2 code, which also 



c/2 . follows chemo-dynamical evolution by accounting for metal enrichment contributed by dif- 

^ i ' ferent stellar populations. We carry out cosmological simulations for a set of galaxy clusters, 

covering the mass range M200 — (0.1 — 2.2) x 10 15 /i -1 A/q. Besides runs not including any 
efficient form of energy feedback, we carry out simulations including three different feedback 
' schemes: ( i) kinetic feedback in the form of galactic winds triggered by supernova explosions; 

(ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (Hi) AGN feed- 
back in which a 'radio mode' is included with an efficient thermal coupling of the extracted 
energy, whenever BHs enter in a quiescent accretion phase. Besides investigating the resulting 
thermal properties of the ICM, we analyse in detail the effect that these feedback models have 
on the ICM metal-enrichment. We find that AGN feedback has the desired effect of quenching 
I star formation in the brightest cluster galaxies at z < 4 and provides correct temperature pro- 

files in the central regions of galaxy groups. However, its effect is not yet sufficient to create 
"cool cores" in massive clusters, while generating an excess of entropy in central regions of 
galaxy groups. As for the pattern of metal distribution, AGN feedback creates a widespread 
enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from 
galactic halos to the inter-galactic medium. This turns into profiles of Iron abundance, Z-p e , 
which are in better agreement with observational results, and into a more pristine enrichment 
of the ICM around and beyond the cluster virial regions. Following the pattern of the relative 
abundances of Silicon and Iron, we conclude that a significant fraction of the ICM enrichment 
is contributed in simulations by a diffuse population of intra-cluster stars. Our simulations 
also predict that profiles of the Zsi/Z-p c abundance ratio do not increase at increasing radii, at 
least out to 0.5R V i r . Our results clearly show that different sources of energy feedback leave 
distinct imprints in the enrichment pattern of the ICM. They further demonstrate that such 
imprints are more evident when looking at external regions, approaching the cluster virial 
boundaries. 

Key words: Cosmology: Theory - Methods: Numerical - X-Rays: Galaxies: Clusters - 
Galaxies: Abundances - Galaxies: Intergalactic Medium 
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1 INTRODUCTION 

High quality data from the current generation of X-ray satel- 
lites (XMM-Newton, Chandra and Suzaku) have now established 
a number of observational facts concerning the thermo-dynamical 
and chemo-dynamical properties of the intra-cluster medium (ICM) 
for statistically representative sets of galaxy clusters: core re- 



gions of relaxed clusters show little eviden ce of gas cooler than 
about a third of the virial temperature (e.g. , IPeterson et alj|200"ll ; 
iBohringer et al.ll2002klSanderson et alj|2006h : temperature profiles 
have negative gradients outside core regions, a trend that ex- 
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2008b); gas entropy is higher than expected from simple self- 
similar scaling properti es of the ICM, n ot only in core regions, but 
also out to -R50CU ( e -g-, 

ISun et alJ 2009. and references therein); ra- 
dial profiles of the Iron abundance show negative gradients, more 
pronounced for relaxed "cool core" clusters, with central values 
of Zyc approaching solar abundance and with a global enrich- 
ment at a level of about 1/3-1/2 Z Fe _ , p (e.g . |Pe Grandi et al 12004 

i. ISnowden et alj \200® 
Werner et al.l 



Vikhlinin et al. 2005; dePlaaetal. 2006; 
Leccardi & Molendil l2008al ; see iMushotzkvl |2004 
20081 for recent reviews). 



Galaxy clusters arise from the collapse of density perturba- 
tions over a scale of ~ 10 h^ 1 comoving Mpc. As such, the 
above observational properties of the ICM come from a non-trivial 
interplay between the underlying cosmological scenario, which 
shapes the large-scale structure of the Universe, and a number of 
astrophysical processes (e.g., star formation, energy and chemi- 
cal feedback from supernovae and AGN) taking place on much 
smaller scales. Cosmological hydrodynamical simulations repre- 
sent the modern instrument with which such complexities can be 
described as the result of hierarchica l assembly of cosmic struc- 
tures (e.g., lBorgani & Kravt sov 2009, for a recent review). Indeed, 
simulations of galaxy clusters reach nowadays a high enough res- 
olution, while including a realistic description of the above men- 
tioned astrophysical processes, for them to provide a coherent 
interpretative framework for X-ray observations. Quite remark- 
ably, simulation predictions for the thermodynamical properties 
of the ICM are in good agreement with observations, at least 
outside the core regions: simulated profiles of gas density and 
temperature mat ch the observed ones at cluster-centric distances 
> O.l flsnn (e.gjLoken et alj|2002l; IBorgani et alj|2004r ~ 



2004 [koncarelli et al. 20061 IPratt et alj|2007l ; iNagai et 



Kay et all 
al .1 120071 : 



Croston et al.ll2008l) : the observed entropy level at -R500 is well re- 



produced by simulations including radiative cooling and star for- 
mation (e.g.. INagai et al.ll2007l ; IBorgani & Vieil lTo09). The situa- 
tion is quite different within cluster cores, where simulations in- 
cluding only stellar feedback generally fail at producing realistic 
cool cores. The two clearest manifestations of this failure are rep- 
resented by the behaviour of temperature and entropy profiles at 
small radii, J; 0.1i?soo- Observations of cool core clusters show 
that temperature profiles smoothly decline toward the centre, reach- 
ing temperatures of about 1/3-1/2 of the maximum value, while the 
entrop y level at the smallest sampled scales is generally very low 
(e.g. Jsun et alj|2009l ; lsanderson et al]|2009ah . On the contrary, ra- 
diative simulations including a variety of models of stellar feedback 
predict steep negative temperature gradients down to the innermost 
resolv ed radii and central entropy levels much higher than observed 
(e.g.. IValdarninil[2003l : iTornatore et ai]|2003l ; IBorgani et ai]|2004 
INagai et alj|20071 ; cf. also lKav et alj|2007h . This failure of simula- 
tions is generally interpreted as due to overcooling, which takes 
place in simulated clusters even when including an efficient su- 
pernova (SN) feedback, and causes an excess of star formation in 
the simulated br ightest cluster galaxies (BCGs; lRomeo et alj2005l ; 
ISaro et al.ll2006l) . 

The generally accepted solution to these shortcomings of sim- 
ulations is represented by AGN feedback. Indeed, the presence of 
cavities in the ICM at the cluster centre is considered as the finger- 



1 Here and in the following we will indicate with the radius encom- 
passing a mean density equal to Ap cr , being p cr = 3H 2 /8irG the critical 
cosmic density. In a similar way, we will indicate with the mass con- 
tained within i?A ■ 



print of the conversion of mechanical energy associated to AGN 
jets into thermal energy (and possibly in a non-thermal content 
of relativistic particles) through shocks (e .g., Mazzotta et al. |200' 
Fabi an et alj2005l;lMcNamara et al.l2 006; Sand ers & FabianhoO' 
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34; 
37; 



McNa mara & Nulsenll2007l for a review). Although analytical 



arguments convincingly show that the energy radiated from gas ac- 
cretion onto a central super-massive black hole (BH) is enough to 
suppress gas cooling, it is all but clear how this energy is ther- 
malised and distributed in the surrounding medium. A likely sce- 
nario is that bubbles of high entropy gas are created at the ter- 
mination of jets. Buoyancy of these bubbles i n the ICM then dis- 
tributes thermal energy ove r larger scales (e.g., |Palla Vecchia et al.l 
|2004 ICattaneo et alj|2007h . Crucial in this process is the stabil- 
ity of the bubbles against gas dynamical instabilities which would 
tend to destroy them quite rapidly. Indeed, detailed numerical sim- 
ulatio ns have demonstra t ed that gas circulation associated to jets 
(e.g. lOmma et al.112004 iBrighenti & Mathews! 120061; iHeinz et al.l 
12006b . gas viscosity, magnetic fiel ds (e.g.,[R uszkowski et aU2007h 
and injection of cosmic rays (e.g.. lRuszkowski et alj|200a) are all 
expected to cooperate in determining the evolution of buoyant bub- 
bles. Although highly instructive, all these simulations have been 
performed for isolated halos and, as such, they describe neither the 
cosmological growth and merging of black holes, nor the hierarchi- 
cal assembly of galaxy clusters. 

ISpringel et all d2005h and |Pi Matteo et al.l J2005h have pre- 
sented a model that follows in a cosmological simulation the evolu- 
tion of the BH population and the effect of energy feedback result- 
ing from gas accretion onto BHs. In this model, BHs are treated as 
sink particles which accrete from the sur rounding gas according to 
a Bondi accretion rate (e.g., Bondi 19 53), with an upper lim it pro- 
vided by the Eddingt on rate (see a lso B ooth ^&" S chavdl2009h . This 
model was used bv lBhattacharva et alj d2008b . who run cosmolog- 
ical simulations to study the effect of AGN feedback on galaxy 
groups. They found that this feedback is effective in reducing star 
formation in central regions and to displace gas towards outer re- 
gions. However, the resulting entropy level within cluster cores is 
higher than observed. In its original version, the energy extracted 
from BH accretion is locally distributed to the gas around the BHs 
according to a SPH-kernel weighting s cheme. This model has been 
modified bv lSiiacki & Sprin gel (2006) who included the possibil- 
ity to inflate high-entropy bubbles in the ICM whenever accretion 
onto the central BH enters in a quiescent "radio mode". The under- 
lying idea of injecting bubbles with this prescription was to pro- 
vide a more realistic description of the effect of jet termination 
on the ICM, although the effect of the jet itself was not included. 
IPuchwein et al. (2008) simulated a set of clusters and groups to 
show that this feedback scheme is able to reproduce the observed 
slope of the relatio n between X-ray luminosity and temperature. 
ISiiacki et alj d2007h showed in a simulation of a single poor galaxy 
cluster that the injection of bubbles is quite effective in suppress- 
ing star formation in central regions. However, also in this case the 
temperature profile in the core regions does not match the observed 
slope, while more realistic temperature profiles can be produced if 
bubble injection is associated t o the injection of a non-thermal pop- 
ulatio n of relativistic particles jPfrommer et alj[2007l ; ISiiacki et al.1 
2008). Although these authors presented results concerning the ef- 
fect of AGN on the ICM thermodynamics, no detailed analysis has 
been so far carried out to study the interplay between AGN and 
chemical enrichment, by including a detailed chemo-dynamical de- 
scription of the ICM. 

Indeed, the amount and distribution of metals in the ICM pro- 
vide an important diagnostic to reconstruct the past history of star 
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formation and the role of gas-dynamical and fee dback processes in 
displacing metals from star forming region s (e.g . jMushotzkvl2004 
iBoi-gani et ai1l2008l : ISchindler & Diaferickooa for reviews). Cos- 
mological chemo-dynamical simulations of galaxy clusters gen- 
erally show that the pre dicted profiles of Iron abundance are 
steeper than observed (e.g..lValdarninill2003l;lTornatore et alj2004l ; 
iRomeo et af]|2005l ; iTornatore et alj|200l iDave et al.ll2008h . with 
an excess of enrichment in the core regions. This is generally inter- 
preted as due to the sa me excess of recent star formation in simu- 
lated BCGs. However, iFabian et alj {2008) showed that an excess 
of recent star formation has also the effect of efficiently locking 
recently-produced metals into stars, thereby preventing a too fast 
increase of the metallicity of the hot diffuse medium. 

All the above analyses are based on different implementa- 
tions of SN energy feedback, while only much less detailed anal- 
yses of ICM metal enrichment have been so far presented by 
also including AGN feedback in co smological si mulations (e.g., 
ISiiacki et ai]|2007l : iMoll et al.|[20071) . For instance, iRoediger etal] 
d2007h showed from simulations of isolated halos that buoyancy of 
bubbles can actually displace a large amount of the central highly 
enriched ICM, thus leading to a radical change of the metallicity 
profiles, or even to a disruption of the metallicity gradients. 

The aim of this paper is to present a detailed analysis of cos- 
mological hydrodynamical simulations of galax y clusters, w hich 
have been carried out with the GADGET-2 code (Springel 2005), by 
combi ning the AGN feedback model described by ISpringel et al.1 
(2005) with the SPH im plementation of chemo-dynamics presented 
bv lTornatore et al.1 J2007h . Besides showing results on the effect of 
combining metallicity-dependent cooling and AGN feedback on 
the ICM thermodynamics, we will focus our discussion on the dif- 
ferent effects that SNe and AGN feedback have on the chemical 
enrichment of the ICM. Although we will shortly discuss the effect 
of different feedback sources on the pattern of star formation in the 
BCG, the results presented in this paper will mainly concern the 
enrichment of the hot diffuse intra-cluster gas. We defer to a forth- 
coming paper a detailed analysis of the effect of AGN feedback on 
the population of cluster galaxies. The scheme of the paper is as 
follows. We present in Section 2 the simulated clusters, and briefly 
describe the implementation of the chemical evolution model in 
the GADGET-2 code along with the SN and AGN feedback mod- 
els. In Section 3 we show our results on the thermal properties of 
the ICM and their comparison with observational results. Section 
4 is devoted to the presentation of the results on the ICM chemi- 
cal enrichment from our simulations and their comparison with the 
most recent observations. We discuss our results and draw our main 
conclusions in Section 5. 



2 THE SIMULATIONS 

2.1 The set of simulated clusters 

Our set of simulated clusters is obtained from nine Lagrangian 
regions extracted from a simulation co ntaining only dark m atter 
(DM) with a box size of 479 h" 1 Mpc jYoshida et aT]|200ll) . per- 
formed for a flat A cold dark matter (CDM) cosmological model 
with relevant parameters Q. m — 0.3, hioo = 0.7, as = 0-9- 
This s et of simulated clusters is described in detail bv lDolag et al.l 
(2008). The regions are centred on nine galaxy clusters: five of them 
have masses M200 — 1-0 x 10 14 h -1 Mq, while the remaining four 
have M200 ^ (1.0 - 2.2) x 10 15 /i _1 M . The regions of three 
largest objects also contain other clusters, so that we have in total 



Table 1. Characteristics of the simulated clusters. Column 1: cluster name. 
Column 2: mass contained within -R200 (units of 10 14 /i —1 Mq). Col- 
umn 3: value of Rqqo (units of h~ x Mpc). Column 4: value of the 
spectroscopic-like temperature within R500, Tsoo- Columns 5 and 6: mass 
of the central BH hosted in the BCG for the AGN1 and AGN2 runs (units of 
10 10 /i _1 M0), respectively. For the g51 cluster, the two additional values 
reported for the mass of the central BH refer to the AGN2(0.8) and AGN2W 
runs (see text) and are indicated with the t and | symbol, respectively. 
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0.64 


1.44 
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0.95 
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0.88 


0.72 


2.03 
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1.13 
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0.83 


0.71 


1.79 


2.54 


0.79 
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0.92 


0.74 


2.03 


3.15 


1.29 
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0.89 


0.73 


2.02 


2.57 


0.99 



Table 2. Resolution of the different runs. Column 2-3: mass of the DM par- 
ticles and initial mass of the gas particles (units of 10* h~ 1 Mq). Column 
4: value of the Plummer-equivalent gravitational force softening at z = 
(units of /i _1 kpc). Column 5: clusters simulated at each resolution. 



hi DM 



Clusters 



Lowres. 11.0 2.03 5.00 gl, g8, g51, g72, g676 
High res. 1.69 0.31 2.75 g676, g914, gl542, g3344, 

g6212 

18 clusters with k'haofL 5 x 10 13 h~ 1 MQ, whose basic character- 
istics are reported in TableQ] 

Each region was re -simulated using the zoomed initial condi- 
tion (ZIC) technique bv lTormen et al.l ( 1997). The relative mass of 
the gas and dark matter (DM) particles within the high-resolution 
region of each simulation is set so that Qbar ~ 0.045 for the den- 
sity parameter contributed by baryons. For each Lagrangian region, 
initial conditions have been generated at a basic resolution and at 
6.5 times higher mass resolution. At this higher resolution, the mass 
of the DM and gas particles are iudm — 1-9 x 10 s h' 1 Mq and 
m gas = 2.8 x 10 7 /i _1 Mq respectively, with Plummer equiva- 
lent softening length for the computation of the gravitational force 
set to £ = 2.75/i _1 kpc in physical units below z = 2 and 
to e — 8.25/i _1 kpc in comoving units at higher redshifts. In 
the lower resolution runs, these values are rescaled according to 

— 1/3 

m DM ■ We decided to simulate the nine Lagrangian regions with 
massive clusters at the lower resolution, while we used the higher 
resolution for the five low-mass clusters. Only the low-mass g676 
cluster was simulated at both resolutions. We provide in Table|2]a 
description of the parameters for the two different resolutions, also 
listing the clusters simulated at each resolution. 



4 Fabjan et al. 



2.2 The simulation code 

Our simulations were performe d using the TreePM-SPH 
GADGET- 2 code dSpringell l2005h . All si mulations include a 
metal licity-dependent radiative cooling dSutherland & Dopita 
1 1993b . h eating from a uniform t ime-dependent ultraviolet back- 
ground jH aardt & Madau 1996) and the effective model by 
ISpringel & Hernquistl 1 2003 ) for the description of star formation. 
In this model, gas particles above a given density are treated as 
multiphase, so as to provide a sub-resolution description of the 
inter-stellar medium. In the following, we assume the density 
threshold for the onset of star formation in multiphase gas particles 
to be riH = 0.1cm -3 in terms of number density of hydrogen 
atoms. Our si mulations also i nclude a detailed model of chemical 
evolution by iTornatore et alj (2007, T07 hereafter). We address 
the reader to T07 for a more detailed description of this model. 
Metals are produced by SNe-II, SNe-Ia and intermediate and 
low-mass stars in the asymptotic giant branch (AGB hereafter). 
We assume SNe-II to arise from stars having mass above 8Mq. 
As for the SNe-Ia, we assume their progenitors to be binary 
systems, whose total mass lies in the range (3-16)M0. Metals 
and energy are released by stars of different mass by properly 
accounting for mass-dependent lifetimes. In this work we as- 
sume the lifetime function proposed by Pad ovani & Matteuccil 
dl993l) . while we assume the standard stellar initial mass function 
(IMF) by ISalpeteri dl955l). We adopt the metallicity-dependent 



stel lar yields by Woosley & Weaverl 1 1995) for SNe-II, the yields 
by Ivan den Hoek & Groenewegenl d 19971) for the AGB and by 
Thieleman n et al] d2003l) for SNe-Ia. The version of the code used 
for the simulations presented here allowed us to follow H, He, C, 
O, Mg, S, Si and Fe. Once produced by a star particle, metals are 
then spread to the surrounding gas particles by using the B-spline 
kernel with weights computed over 64 neighbours and taken to be 
proportional to the volume of each particle (see T07 for detailed 
tests on the effect of choosing different weighting schemes to 
spread metals). 



2.3 Feedback models 

In the simulations presented in this paper we model two different 
sources of energy feed back. The first one is the kin etic feedback 
model implemented by Springel & Hernquist] d2003h . in which en- 
ergy released by SN-II triggers galactic winds, with mass up- 
load rate assumed to be proportional to the star formation rate, 
Mw = riMi,. Therefore, fixing the parameter 77 and the wind ve- 
locity vw amounts to fix the total energy carried by the winds. 
In the following, we assume r\ = 2 for the mass-upload parameter 
and vw = 500 km s _1 for the wind velocity. If each SN-II releases 
10 51 ergs, the above choice of parameters corresponds to assuming 
that SNe-II power g alactic outflows with nearly unity efficiency for 
a Salpeter IMF (see lSpringel & Hernauisll2003h . 

Furthermore, we include in our simulations the effect of 
feedback energy released by gas accretion onto super-massive 
black holes (B Hs), following the scheme original ly introduced by 
Springel et akl (2005, SDH05 here after; see also |Pi Matteo et alj 
2005Ll2008hlBooth & Schavel2009h . We refer to SDH05 for a more 
detailed description. In this model, BHs are represented by colli- 
sionless sink particles of initially very small mass, that are allowed 
to subsequently grow via gas accretion and through mergers with 
other BHs during close encounters. During the growth of structures, 
we seed every new dark matter halo above a certain mass threshold 
Mth, identified by a run-time friends-of-friends algorithm, with a 



central BH of mass 10 5 h' 1 Mq, provided the halo does not con- 
tain any BH yet. For the runs at the lower resolution the value of 
the halo mass assumed to seed BHs is M t h = 5 x 10 10 h' 1 Mq, 
so that it is resolved with about 40 DM particles. At the higher 
resolution the halo mass threshold for BH seeding decreases to 
M t h = 10 10 hT 1 Mq, so as to resolve it with approximately the 
same number of particles. We verified that using in the higher- 
resolution runs the same value of Mth as in the low-resolution runs 
causes the effect of BH accretion to be shifted toward lower red- 
shift, since its onset has to await the formation of more massive 
halos, while leaving the final properties of galaxy clusters almost 
unaffected. 

Once seeded, each BH can then grow by local gas accretion, 
with a rate given by 



Mbh = min (Mb, Me 



(1) 



or by merging with other BHs. Here M b is the accretion rate esti- 
mated with the Bondi-Hoyle-Lyttleton for mula dHovle & Lyttletor] 
ll939l;lBondi & Hovldfl944llBondilll952|) . while M Edd is the Ed- 
dington rate. The latter is inversely proportional to the radiative 
efficiency e,-, which gives the radiated energy in units of the energy 
associated to the accreted mass: 



(2) 



M BH c 2 

Following Sprin gel et all {2005), we use e r = 0.1 as a reference 
value, which is typi cal for a radiatively efficien t accretion onto a 
Schwartzschild BH ( Shak ura & Svunvaevlll973h . The model then 
assumes that a fraction e/ of the radiated energy is thermally 
coupled to the surrounding gas, so that Ef cc d = e r e/Msirc 2 
is the rate of provided energy feedback. In stan dard AGN feed- 
back i mplementation we use e/ = 0.05 following |Pi Matteo et alj 
(2005), who were able with this value to reproduce the observed 
Mbh — a relation b etween bulge velocity dispersion and mass of 
the hosted BH (e.g.. lMagorrian et al.l l 1998). This choice was also 
found to be consistent with the value required in semi-analytical 
models to explain the evolution of the number density of quasars 
dWvithe & Loebll2003l) . 

Gas swallowed by the BH is implemented in a stochastic way, 
by assigning to each neighbour gas particle a probability of con- 
tributing to the accretion, which is proportional to the SPH kernel 
weight computed at the particle position. Differently from SDH05, 
we assume that each selected gas particle contributes to the accre- 
tion with 1/3 of its mass, instead of being completely swallowed. 
In this way, a larger number of particles contribute to the accretion, 
which is then followed in a more continuous way. We remind that 
in the SDH05 scheme, this stochastic accretion is used only to in- 
crease the dynamic mass of the BHs, while their mass entering in 
the computation of the accretion rate is followed in a continuous 
way, by using the analytic expression for Mbh- Once the amount 
of energy to be thermalised is computed for each BH at a given 
time-step, one has to distribute this energy to the surrounding gas 
particles. In their original formulation, SDH05 distributed this en- 
ergy using the SPH kernel. 

Besides following this standard implementation of the AGN 
feedback, we also follow an alternative prescription, which differ 
from the original one i n two aspects . 

Firstly, following Siiac ki et alj d2007l) . we assume that a tran- 
sition from a "quasar" phase to "radio" mode of the BH feedback 
takes place whenever the accretion rate falls below a given limit 
(e.g.. lChurazov et ai]|2005L and references therein), corresponding 
to Mbh / MEdd < 10 -2 , which implies an increase of the feed- 
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back efficiency to e/ = 0.2. At high redshift BHs are characterised 
by high accretion rates and power very luminous quasars, with only 
a small fraction of the radiated energy being thermally coupled to 
the surrounding gas. On the other hand, BHs hosted within very 
massive halos at lower redshift are expected to accrete at a rate 
well below the Eddington limit, while the energy is mostly released 
in a kinetic form, eventually thermalised in the surrounding gas 
through shocks. Secondly, instead of distributing the energy using 
a SPH kernel, we distribute it using a top-hat kernel, having radius 
given by the SPH smoothing length. In order to avoid spreading the 
energy within a too small sphere, we assume a minimum spread- 
ing length of 2/i _1 kpc. The rationale behind the choice of the 
top-hat kernel is to provide a more uniform distribution of energy, 
thus mimicking the effect of inflating bubbles in correspondence 
of the termination of the AGN jets. It is well known that a num- 
ber of physical processes need to be adequately included for a fully 
self-consistent description of bubble injection and buoyancy: gas- 
dynamical effects related to jets, magnetic fields, viscosity, thermal 
conduction, injection of relativistic particles. 

In particular, a number of studies based on simulations of iso- 
lated halos (e.g., IOmmaetalJl2004 iBrighenti & Mathews! 1 2006) 
have pointed out that gas circulation generated by jets provides an 
impor tant contribution for the stabilization of cooling flows (see 
also lHeinz et alj2006l , for a cosmological simulation of cluster for- 
mation including jets). In its current implementation, the model of 
BH feedback included in our simulations neglects any kinetic feed- 
back associated to jets. Based on an analytical model, Pope ( 2009) 
computed the typical scale of transition from kinetic to thermal 
feedback regime for AGN in elliptical galaxies and clusters. As 
a result, he found that the effect of momentum carried by jets can 
be neglected on scales <^ 20 kpc, the exact value depending on the 
local conditions of the gas and on the injection rate of kinetic and 
thermal energy. In order to compare such a scale to that actually 
resolved in our simulations, we remind the reader that SPH hydro- 
dynamics is numerically converged on scales about 6 times larger 
than the Plummer-equiva lent softening scale for gravitational force 
(e.g., lBorganie7al]|2002l) . Owing to the values of the gravitational 
softening reported in Table 2, the scales resolved in our simulations 
are not in the regime where kinetic feedback is expected to domi- 
nate, thus justifying the adoption of a purely thermal feedback. As 
a further test, we have computed the radius of the top-hat kernel 
within which energy is distributed around the central BHs in our 
simulated clusters. As a few examples, we found at z — this ra- 
dius to be 21 kpc and 23 kpc for the AGN2 runs of the g72 and g676 
clusters, respectively. This implies that we are in fact distributing 
thermal energy over scales where kinetic feedback should not be 
dominant. 

In view of the difficulty of self-consistently including the co- 
operative effect of all the physical processes we listed above, we 
prefer here to follow a rather simplified approach and see to what 
extent [results on] the final results of our simulations are sensitive 
to variations in the implementation of the BH feedback model. 

In summary, we performed four series of runs, corresponding 
to as many prescription for energy feedback. 

(a) No feedback (NF hereafter): neither galactic winds nor AGN 
feedback is included. 

( b) Galactic winds (W hereaft er) included by following the model 
by Springel & Hernquist] d2003l) . with v w — 500 km s~ 1 and r\ = 2 
for the wind mass upload. 

(c) Standard implementation of AGN feedback from BH accre- 
tion (AGN1 hereafter), using e/ = 0.05 for the feedback efficiency. 



(d) Modified version of the AGN feedback (AGN2 hereafter), 
with feedback efficiency increasing from ey = 0.05 to e/ = 0.2 
when Mbh /Meu < 10~ 2 , and distribution of energy around the 
BH with a top-hat kernel. 

In order to further explore the parameter space of the consid- 
ered feedback models we also carried out one simulation of the 
g5 1 cluster based on the AGN2 scheme, but with the feedback effi- 
ciency increased to €j = 0.8 (AGN2(0.8) hereafter). We also note 
that our simulations include either winds triggered by SN explo- 
sions or AGN feedback. While this choice is done with the purpose 
of separate the effects of these two feedback sources, we expect in 
realistic cases that both AGN and SN feedback should cooperate 
in determining the star formation history of galaxies. In order to 
verify the effect of combining the two feedback sources, we car- 
ried out one simulation of the AGN2 scheme for the g51 cluster, 
in which also galactic winds with a velocity v w = 300 kms -1 are 
included (AGN2W hereafter). 

Before starting the presentation of the results on the thermal 
and enrichment properties of the ICM, we briefly comment on the 
results concerning the mass of the central black holes in the simu- 
lations including AGN feedback, and the star formation rate (SFR) 
history of the brightest cluster galaxies (BCGs). 

Looking at Table 1, we note that our simulations predict 
rather large masses for the super-massive BHs hosted in the central 
galaxies. Quite interestingly, the AGN2 runs generally produce BH 
masses which are smaller, by a factor 3-5, than for the AGN1 runs. 
This demonstrates that including the more efficient "radio mode" 
for the feedback and distributing the energy in a more uniform 
way has a significant effect in regulating gas accretion. Although 
BCGs are known to host BHs which are more massive than ex- 
pected for norma l early type galaxies of comparable mass (e.g., 
lLauer et alj|2007h . the BH masses from our simulations seem ex- 
ceedingly large, also for the AGN2 model. For instance, the BH 
mass hosted by M87, within the re latively poor Virgo c luster, is 
m B H ^ (3 - 4) x 1O 9 M (e.g., iRaffertv et alj|2006r) . This is 
about a factor 3-5 smaller than found in our AGN2 runs for clus- 
ters of comparable mass, M200 ^2x 10 14 /i _1 M . We note that 
increasing the radio-mode feedback efficiency from ef = 0.2 to 0.8 
(AGN2(0.8) run) reduces the mass of the central black hole in the 
g51 cluster from m BH ^ 16.5 x 1O 9 M to ~ 6.2 x 1O 9 M (see 
Table 1). This suggests that a highly efficient thermalization of the 
energy extracted from the BH is required to regulate gas accretion 
to the observed level. Quite interestingly, we also note that adding 
the effect of galactic winds in the AGN2 scheme (AGN2W run) 
only provides a marginal reduction of the final mass of the central 
BH. Therefore, although galactic winds can play a significant role 
in regulating star formation within relatively small galaxies, they 
are not efficient in decreasing gas density around the largest BHs, 

so as to suppress their accretion rate. 

As for the comparison with previous analyses, Siiac kTet al.l 

(2007) performed a simulation of the same g676 cluster included 
in our simulation set at the lower resolution, for their feedback 
scheme based on the injection of AGN driven bubbles. They found 
that the final mass of the central BH is tubh — 6 x 10 9 /i -1 M . 
This value is about 30 per cent lower than what we find for the 
AGN2 ru ns of g676. In ord er to compare more closely with the 
result by ISiiacki et alj d2007t) . we repeated the run of the g676 
cluster by switching off the metallicity-dependence of the cool- 
ing function. As a result, BH accretion proceeds in a less effi- 
cient way, and the resulting central BH mass in this case drops to 
~ 3.5 x W 9 h~ 1 M e . 
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As for the SFR history of the BCG, this is estimated by identi- 
fying first all the stars belonging to the BCG at z = 0. This has bee n 
done by running the SKID group-finding algorithm (Stadel 2001), 
using the same procedure described bv lSaro et al.l (2006). After tag- 
ging all the star particles belonging to the BCG, we reconstruct the 
SFR history by using the information on the redshift at which each 
star particle has been spawned by a parent gas particle, according 
to the stochastic algorithm of star forma t ion im plemented in the ef- 
fective model by Spri ngel & Hernquistl ( 120030 . The resulting SFR 
histories for the BCG of the g5 1 clusters are shown in the left panel 
of Figure Q] 

As expected, the run with no efficient feedback (NF) produces 
the highest star formation at all epochs: SFR has a peak at z ~ 4 
and then drops as a consequence of the exhaustion of gas with short 
cooling time. Despite this fast reduction of the SFR, its level at 
z = is still rather large, ~ 400 h -1 M@ yr" 1 . As for the run 
with winds, it has a reduced SFR since the very beginning, owing 
to the efficiency of this feedback scheme in suppressing star for- 
mation within small galaxies which form already at high redshift. 
Also in this case, the peak of star formation takes place at z ~ 4, 
but with an amplitude which is about 40 per cent lower than for the 
NF runs and with a more gentle decline afterwards. As for the runs 
with AGN feedback, their SFR history is quite similar to that of the 
W run down to z ~ 5. Star formation is then suddenly quenched 
at z^S 4. We note in general that SFR for the AGN1 model lies 
slightly below that of the AGN2 model, as a consequence of both 
the different way of distributing energy and, at low redshift, of the 
inclusion of the radio mode assumed in the quiescent BH accretion 
phase. Suppression of the SF at relatively low redshift is exactly 
the welcome effect of AGN feedback. At z — 0, the resulting SFR 
is of about 7OM0yr _1 for both models, a value which is closer to, 
although still higher than the typical values o f SFR observed in th e 
BCGs of clusters of comparable mass. (e.g.. lRaffertv"eUd1 l2006). 
Quite interestingly, increasing the feedback efficiency to e/ = 0.8 
in the AGN2(0.8) run does not significantly affect the level of low- 
redshift star formation, while it suppresses star formation by only 
~ 10 per cent around the peak of efficiency. Therefore, while a 
higher efficiency is indeed effective in suppressing gas accretion 
onto the central BH, a further reduction of star formation associ- 
ated to the BCG should require a different way of thermalizing the 
radiated energy. As for the AGN2W run, its SFR at high redshift is 
lower than for the AGN2 run, due to the effect of winds, while no 
significant change is observed at z^S 3. 

In the right panel of Fig[T|we plot the stellar mass found in the 
BCG at z = that is formed before a given redshift. According 
to this definition, this quantity is the integral of the SFR plotted in 
the right panel, computed from a given redshift z to infinity. This 
plot clearly shows the different effect that winds and AGN feed- 
back have in making the BGC stellar population older. In the NF 
and W runs, the redshift at which 50 per cent of the BCG stellar 
mass was already in place is z$o — 2. This indicates that even an 
efficient SN feedback is not able to make the stellar population of 
the BCG older. On the contrary, the effect of AGN feedback takes 
place mostly at relatively low reshift. As a consequence the age of 
the BGC stellar population increases, with 250 — 3.0 — 3.6, almost 
independent of the detail of the AGN feedback scheme. 



3 THERMAL PROPERTIES OF THE ICM 

Galaxy clusters are identified in each simulation box by running 
first a friends-of-friends (FOF) algorithm over the high resolution 



DM particles, using a linking length of 0.15 in terms of the mean 
inter-particle separation. Within each FOF group, we identify the 
DM particle having the minimum value of the gravitational poten- 
tial and take its position to correspond to the centre of the cluster. 
All profiles are then computed starting from this centre. The small- 
est radius that we use to compute profiles encompasses a minimum 
number of 100 SPH particles, a criterion that gives numerically 
converged results for profiles of gas density and t emperature in 
non-radiative simulations (e.g.. iBorgani et al]l2002h . As we shall 
see in the following, profiles computed with this criterion will ex- 
tend to smaller radii for those runs which have a higher gas density 
at the centre, while stopping at relatively larger radii for the runs in- 
cluding AGN feedback, which are characterised by a lower central 
gas density. 

3.1 The luminosity-temperature relation 

The relation between bolometric X-ray luminosity, Lx and ICM 
temperature, T, provided one of the first evidences that non- 
gravitational effects determ ine the thermo-dynamical properties of 
the ICM (e.g.. I Voi3 12005L for a review). Its slope at the scale of 
clusters is observed to be Lx oc T a with a ~ 2.5-3 (e.g., [Hornet] 
120011 : IPratt et"alll200gt) . and possibly even ste eper or with a larger 
scatte r at the scale of galaxy groups (e.g., lOsmond & PonmarJ 
2004). These results are at variance with respect to the prediction, 
q = 2, of self-similar mo dels based o nly on the effect of gravita- 
tional gas accretion (e.g., Kaiser 1986). Attempts to reproduce the 
observed Lx-T relation with hydrodynamic simulations of clus- 
ters have been pursued by several groups (see Borga ni et al. I l2008l. 
for a recent review). Simulations of galaxy clusters including the ef- 
fect of star formation and SN feedback in the form of energy-driven 
galactic winds produce results which are close to observations 
at the scale of clus ters, while general ly producing too luminous 
galaxy groups (e.g.. lBorgani et af]|2004l) . Although a closer agree- 
ment with observations at the scale of groups c an be obtained b y 
using SN-triggered momentum-driven winds l lDave et all 12008). 
there is a general consensus that stellar feedback can not reproduce 
at the same time both the observed Lx-T relation and the low star 
forma tion rate observed in central cluster galaxies. IPuchwein et al.l 
d2008l) presented results on the Lx-T relation for simulations of 
galaxy clusters which included the bubble -driven AGN feedback 
scheme introduced by ISiiacki et al.l J2007t) . They concluded that, 
while simulations not including any efficient feedback (in fact, 
quite similar to our NF runs) produce overluminous objects, their 
mechanism for AGN feedback is efficient in suppressing the X-ray 
luminosity of clusters and groups at the observed level. 

We present here our results on the Lx-T r elation, keeping 
in min d that our simulations differ from those by IPuchwein et al] 
(2008) both in the details of the implementation of the AGN feed- 
back scheme and in the treatment of the metallicity dependence of 
the cooling function. In the left panel of Figure [2] we show the re- 
sults for our runs based on SN galactic winds (W runs) and for the 
runs not including any efficient feedback (NF runs). Filled sym- 
bols refer to the main halo of each simulated Lagrangian region, 
while open circles are for the "satellites". Although we find sev- 
eral satellites having a temperature comparable to those of the low- 
mass main halos (see also Table 1), we remind that these satellites 
are described with a mass resolution six times lower than that of the 
low-mass main halos. We note that the NF runs provide a Lx-T re- 
lation which is not far from the observed one, especially at the scale 
of group s. The reason for this cl oser agreement, with respect to the 
result by IPuchwein et alj d2008l) lies in the fact that these authors 
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Figure 1. Left panel: the history of star formation rate of the brightest cluster galaxy (BCG) in the g51 cluster for the runs with no feedback (NF, blue short- 
dashed), with galactic winds (W, green dashed), with standard AGN feedback (AGN1, red dot-dashed) and with modified AGN feedback (AGN2, light-blue 
solid). Also shown are the results for the AGN2 run with higher radio-mode BH feedback efficiency (AGN2(0.8), purple long-dashed) and for the AGN2 run 
also including galactic winds (AGN2W, magenta dotted). Right panel: the stellar mass found in the BCG at z = that was formed before a given redshift. 



adopted a cooling function computed for zero metallicity. Includ- 
ing the contribution of metal lines to the radiative losses increases 
cooling efficiency and, therefore, gas removal from the hot X-ray 
emitting phase. However, the price to pay for this reduction of X- 
ray luminosity is that a far too large baryon fraction is converted 
into stars within clusters (see below). Quite paradoxically, we also 
note than including an efficient feedback in the form of galactic 
winds turns into an increase of X-ray luminosity. This is due to the 
fact that this feedback prevents a substantial amount of gas from 
cooling, without displacing it from the central cluster regions, thus 
increasing the amount of X-ray emitting ICM. 

In order to asses the effect of the different resolution used for 
large and small clusters, we carried out runs of g676 at the same 
lower resolution of the massive clusters. The results are shown in 
Fig. [2] with black diamonds connected with arrows to the corre- 
sponding higher resolution result. We note that resolution effects 
go in opposite directions for the NF and W runs. In fact, in the 
absence of winds, the runaway of cooling with resolution removes 
from the hot phase a larger amount of gas, thus decreasing X-ray 
luminosity. On the contrary, higher resolution allows a more ac- 
curate description of kinetic feedback, which starts heating gas at 
higher redshift. As a consequence, radiative losses are compensated 
at higher resolution for a larger amount of diffuse baryons, which 
remain in the hot phase, thereby increasing the X-ray luminosity. 

As for the runs with AGN feedback, it is quite efficient in 
decreasing X-ray luminosity at the scale of galaxy groups, thus 
well recovering the observed Lx-T relation. This conclusion holds 
for both implementations of the AGN feedback, which have rather 
small differences. Therefore, although the AGN2 scheme is more 
efficient in regulating the growth of the central BHs hosted within 
the BCGs, it has a marginal effect on the X-ray luminosity. These 
results are in line with those found bv lPuchwein et al] J2008h . who 
however used the injection of heated bubbles to distribute the en- 
ergy extracted from the BH accretion. This witnesses that the feed- 
back energy associated to gas accretion onto super-massive BHs is 



indeed able to produce a realistic Lx-T relation, almost indepen- 
dent of the detail of how the energy is thermalised in the surround- 
ing medium. 

3.2 The entropy of the ICM 

Entropy level in central regions of clusters and groups is consid- 
ered another fingerprint of the mechanisms which determine the 
thermodynamical history of the ICM. Early observational r esults 
on the presence of entropy cores (e.g.. IPonman et aU 1 19991) have 
been more recently r evised, in the light o f higher quality data 
from Chandra (e.g.. ICavagnolo et al.l l2009h and XMM-Newton 
(e.g. Jjohnson et alj|2009l) observations. These more recent analy- 
ses show that entropy level of clusters and groups at R500 is higher 
than predicted by non-radiative simulations, with entropy profiles 
for relaxed systems co ntinuously decre asing down to the smallest 
resolved radii (see also lSun et al.l l2009). 

While radiative simulations of galaxy clusters have been gen- 
erally shown to reproduce obse rvationalr results at -R500 (e.g., 
iNagai et alj2007l ;l Dave et al . 2008), they g enerally predict too low 
entropy levels at smaller radii. For instance. lBorgani & Viell ( 2009) 
have shown that the entropy at -R2500 can be reproduced by re- 
sorting to a fairly strong pre-heating at z — 4. However, this pre- 
heating must target only relatively overdense regions, to prevent the 
creation of too la rge voids in the s tructure of the Lyman-a forest at 
z ~ 2 (see also lShang et~ail 20071) . 

In the following we use the standard definition of entropy, 
which is usuall y adopted in X-ray studies of galaxy clusters (e.g., 
ISun et alj|2009l) : 

K A = 4^, (3) 

where Ta and n e ,A are the values of gas temperature and electron 
number density computed at Ra. As for the temperature, it is com- 
puted by following the prescription of spectroscopic like temper- 
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T S| «R 500 ) [keV] '° ' T sl «R 500 ) [keV] 

Figure 2. Relation between X— ray luminosity and temperat ure for simulated (col oured symbols) a nd observed (grey points with err orbars) clusters and groups . 
Observational data points are from lArnaud & Evrard fl999» (grey diamonds) and lPratt et ai]d2009l) (grey squares) for clusters, from Osmond & Ponma 
(grey circles) for groups and from Horner 1200 ll) (grey triangles). Data from simulations were computed inside Rsoo- Left panel: results for the no feedback 
case (NF, blue triangles) and for the case with galactic winds (W, green squares). Right panel: results for the runs with standard AGN feedback (AGN1, red 
circles) and with the modified AGN feedback scheme, where also a radio-mode regime is included (see text, AGN2; cyan circles). For each series of runs, 
filled and open symbols refer to the main halo within each resimulated Lagrangian region and to "satellite" halos respectively. Black diamonds refer to the 
runs at 1 X resolution, with arrows pointing to results at 6 X higher resolution for the g676 cluster. 



ature introduced bv lMazzotta et alj J2004I) . This definition of tem- 
perature has been shown to accurately reproduce, within few per- 
cents, the actual spectroscopic temperature obtained by fitting spec- 
tra of simulated clusters with a single-temperature plasma model, 
within the typical energy bands where detectors on-board of present 
X-ray satellites are typically sensitive. 

We show in Figure[3]the com parison between our simulations 
and observational data on groups dSun et al ] |2009h and on clusters 
( VikM inin et al.lr2009h for the relation between entropy and tem- 
perature at i?5oo and R2500 (upper and lower pa nels, respectively) . 
In order to reproduce the procedure adopted bv lSun etaf. (2009), 
we compute the spectroscopic-like temperature of simulated clus- 
ters by excluding the core regions within 0.157?5oo- As for the runs 
with no efficient feedback (NF) we note that they produce entropy 
levels, at both -R500 and -R2500, which are close to the observed 
ones. This result can be explained in the same way as that found 
for the Lx-T relation: overcooling, not balanced by an efficient 
feedback mechanism, removes a large amount of gas from the X- 
ray emitting phase, while leaving in this hot phase only relatively 
high entropy gas, which flows in from larger radii as a consequence 
of lack of central pressure support. 

Including winds (W runs) has the effect of increasing the 
amount of low-entropy gas, which is now allowed to remain in the 
hot phase despite its formally short cooling time, thanks to the con- 
tinuous heating provided by winds. As a result, entropy decreases 
at both radii, for the same reason for which X-ray luminosity in- 
creases. As for the runs with AGN feedback, its effect is almost 
negligible for massive clusters. On the other hand, AGN feedback 
provides a significant increase of the entropy level in poor systems, 
an effect which is larger at smaller cluster-centric radii. The result- 
ing entropy is higher than indicated by observational data. Together 



with the result on the Lx-T relation, this result shows that a rather 
tuned energy injection is required, which must be able to suppress 
X-ray luminosity by decreasing gas density, while at the same time 
reproducing the low entropy level measured at small radii. 

3.3 Temperature profiles 

A number of comparisons between observed and simulated temper- 
ature profiles of galaxy clusters have clearly demomnstrated that a 
remarkable agreement exists at relatively large radii, 0.2i?i8o, 
where the effect of cooling is relatively unimportant. While this re- 
sult, which holds almost indepe ndently of the physical processes 
included in the simulations (e.g., [Token et a l. 2002; B organi et al.1 
|2004 iKav et al.ll2007l ; iPratt et alj|2007l : lNagai etal]|2007t) . should 
be regarded as a success of cosmological simulations of galaxy 
clusters, the same simulations have much hard er time to predict 
realistic profiles within cool-core regions (e.g jBorgani et alj2008l 
for a recent review). In this regime, radiative simulations system- 
atically produce steep negative temperature profiles, at variance 
with observations, as a consequence of the lack of pressure sup- 
port caused by overcooling. This further demonstrates that a suit- 
able feedback mechanisms is required to pressurise the gas, so as 
to prevent overcooling and turning the temperature gradients from 
negative to positive in the core regions. While feedback associated 
to SNe has been proved not to be successful, AGN feedback is gen- 
erally considered as a likely solution for simulations to produce re- 
alistic cool cores. Based on simulations of one relatively low-mass 
cluster. [siiacki et all J2OO8) found that AGN feedback can provide 
reasonable temperature profiles only if a population of relativistic 
particles is injected along with thermal energy in inflated bubbles. 
We present in Figure|4]the comparison between simulated and 
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Figure 3. Relation between entropy and tempe rature for our sim ulated clusters (coloured circ les and squares) and for the observational data points at -R500 
(upper panels) and -R2500 (lower panels) from lSun et alj J2009I) (black filled diamonds) and I Vikhlinin et all 120091) (black open diamonds). Left and right 
panels show results for the nine central clusters for the runs without (NF: blue squares; W: green circles) and with AGN feedback (AGN1: red squares; AGN2: 
cyan circles), respectively. For a fair comparison with observations, spectroscopic-like temperatures of the simulated clusters are computed by excluding the 
regions within 0.15iJsoo. 



observed temperature profiles for galaxy clusters with T^t 3 keV 
(left panel) and for poorer clusters and groups with T 1 ^ 3 keV (right 
panel). Obse rvational results are taken from iLeccardi & Molendil 
(2008b) and lSunetalJ d2009h for rich and poor systems, respec- 
tively. As for rich clusters, none of the implemented feedback 
scheme is capable to prevent the temperature spike at small radii, 
while all models provide a temperature profile quite similar to the 
observed one at B~Sl 0.3i?igo- The situation is different for groups. 
In this case, both schemes of AGN feedback provide results which 
go in the right direction. While galactic winds are not able to signif- 
icantly change the steep negative temperature gradients, the AGN1 



and the AGN2 feedback schemes pressurise the ICM in the central 
regions, thus preventing adiabatic compression in inflowing gas. 
From the one hand, this result confirms that a feedback scheme 
not related to star formation goes indeed in the right direction of 
regulating the thermal properties of the ICM in cool core regions. 
On the other hand, it also demonstrates that the schemes of AGN 
feedback implemented in our simulations do an excellent job at the 
scale of galaxy groups, while they are not efficient enough at the 
scale of massive clusters. 
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Figure 4. Comparison between the temperature profiles for simulated and observed clusters with 3 keV (left panel) and for groups with T< 3 keV (right 
panel). In each panel, different lines corresponds to the average simulated profiles computed over the four main massive clusters in the left panel and over 
the five main low-mass clusters in the right panel, for the different sets of runs: no feedback (NF, blue short dashed), galactic winds (W, green long dashed), 
standard AGN feedback (AGN1, red dot-dashed), modified AGN feedback (AGN2, cyan solid). For reasons of clarity, we show with error bars the r.m.s. scatter 
over the ensemble of simulated clusters only for the AGN 1 runs. Observ ational data points for clusters in the left panel are taken from Leccardi & Molendi 
while those for groups in the right panel are from lSun et al] <2009l) . 
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Figure 5. C omparison of the gas f raction within -R500 in simulations (coloured circles) and observational data from Chandra data (diamonds with errorbars) 
analysed bv lVikhlinin et alj 120061) . Left panel: results for the no feedback (NF) runs (dark blue) and for the runs with galactic winds (W, light green). Right 
panel: results for the runs with the standard AGN feedback (AGN1, dark red) and with modified AGN feedback (AGN2, light cyan). The horizontal dotted line 
marks the cosmic baryon fraction assumed in the simulations. 
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Figure 6. Comparison of the star fraction within -R500 in simulations (coloured circles) and observations (squares with errorbars). Left panel: results for the 
no feedback (NF) runs (dark blue) and for the runs with galactic winds (W, light green). Right panel: results for the runs w ith the standard AGN feedback 
(AGN1, dark red) and with modified AGN feedback (AGN2, light cyan). Observational points are from Gonz alez et alj 120071) where stellar mass includes the 
brightest cluster galaxy (BCG), intra-cluster light (ICL) and galaxies within -R500 ■ 



3.4 The gas and star mass fractions 

The inventory of baryons within galaxy clusters represents an im- 
portant test to understand both the efficiency of star formation and 
how the gas content is affected by feedback mechanisms. In gen- 
eral, the difficulty of regulating gas cooling i n cluster simulations 
causes a too large stellar mass fraction (e.g., Borgani et al]|2004l; 



iKav et al.ll2007l : iNagai et alj2007l: iDave et al.ll2008l) . which in turn 
should correspond to a too low fraction of gas in the diffuse ICM. 

We show in Figure [5] the comparison between simulation re- 
sults and observational data on the mass fraction o f hot gas as a 
function of temperature from lVikhlinin et al, I d2006h . Gas in simu- 
lated clusters is assigned to the hot phase if it is not associated to 
multi-phase gas particles and if its temperature exceeds 3 x 10 4 K. 
A comparison between the runs with no feedback (NF) and with 
galactic winds (W) shows that the latter are characterised in general 
by larger f gas values. This is in line with the results on the Lx-T 
relation and confirms that winds are effective in suppressing cool- 
ing, thus increasing the hot baryon fraction. While there is a rea- 
sonable agreement at the scale of poor clusters, simulations show a 
weak trend with temperature, with f gas for the hotter systems hav- 
ing values well below the observed ones. As for AGN feedback, 
it has the effect of increasing f gaa for the most massive clusters, 
although the resulting gas fraction is still below the observational 
level by about 30 per cent. On the contrary, at the scale of poor clus- 
ters the effect of AGN feedback is that of decreasing f gas below the 
observational limit. Indeed, while in rich systems the effect of AGN 
feedback is that of reducing overcooling, thereby leaving a larger 
amount of gas in the hot phase, in poor systems it is so efficient as 
to displace a large amount of gas outside the cluster potential wells. 

Our result on the low value of f gas at the scale of rich clusters, 
even in the presence of AGN feedback, is in line with the some- 
what low value of Lx seen in Fig. [2] However, thi s resul t is in 
disagreement with that presented by Puchwein et alj | |2008|) . who 



showed instead a good agreement at all temperatures between their 
simulations including AGN feedback and observational data. There 
may be two reasons for this difference. First ly, the scheme to injec t 
AGN-driven high-entropy bubbles used bv lPuchwein et al] (2008) 
could provide a more efficient means of stopping cooling in central 
cluster regions, at the same time preve nting excessive gas rem oval 
in low-mass systems. Secondly, unlike IPuchwein et al.1 ( 2008) wc 
include the dependence of metallicity in the cooling function. As al- 
ready discussed, this significantly enhances cooling efficiency and, 
therefore, the removal of gas from the hot phase. In order to verify 
the impact of this effect, we repeated the AGN2 run of g51, by as- 
suming zero metallicity in the computation of the cooling function. 
As a result, we find that f gas increases from 0.09 to 010. From the 
one hand, this result implies that the more efficient gas accretion 
onto BHs, due to metal-cooling, provides a stronger energy feed- 
back which, in turn, balances the higher cooling efficiency. On the 
other hand, it also implies that the m ain reason for the difference 
with respect to lPuchwein et al] J2OO8) should be rather ascribed to 
the different way in which energy associated to BH accretion is 
thermalised in the surrounding medium. 

As for the behaviour of the mass fraction in stars, we show in 
Figure[6]the comparison between o ur simulations and observational 
results from [Gonzalez et al] d2007t) , who also included in the stel- 
lar budget the contr ibution from diffuse intra-cluster stars (see also 
iGiodini et alj2009r) . These results confirm that none of our simula- 
tions are able to reproduce the observed decrease of f s tar with in- 
creasing temperature. Overcooling is indeed partially prevented in 
the presence of winds and, even more, with AGN feedback. How- 
ever, while simulation results for poor clusters are rather close to 
observations, overcooling in massive clusters is only partially al- 
leviated by AGN feedback, with values of fstar which are larger 
than the observed ones by a factor 2-3. 

In summary, the results presented in this section demonstrated 
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that AGN feedback has indeed a significant effect in bringing the 
Lx-T relation and the entropy level of the ICM closer to observa- 
tional results, while regulating cooling in the central regions. How- 
ever, the effect is not yet large enough to produce the correct tem- 
perature structure in the cool core of massive clusters and, corre- 
spondingly, the correct share of baryons between the stellar and the 
hot gas phase. 



4 METAL ENRICHMENT OF THE ICM 

The X-ray spectroscopic studies of the content and distribution of 
metals in the intra-cluster plasma provides important information 
on the connection between the process of star formation, taking 
place on small scales within galaxies, and the processes which de- 
termine the thermal properties of the ICM. The former affects the 
quantity of metals that are produced by different stellar popula- 
tions, while the latter gives us insights on gas-dynamical processes, 
related both to the gravitational assembly of clusters and to the 
feedback mechanisms that displace metal-enriched gas from star 
forming regions. 

The detailed model of chemical evolution included in the 
GADGET-2 code by T07 allows us to follow the production of 
heavy elements and to study how their distribution is affected by 
the adopted feedback schemes. The analysis presented in this sec- 
tion is aimed at quantifying the different effects that galactic out- 
flows triggered by SN explosions and AGN feedback have on the 
enrichment pattern of the ICM. We present results on the Fe dis- 
tribution and the corresponding abundance profiles, the enrichment 
age within clusters and groups, the relation between global metal- 
licity and ICM temperature, and the relative abundance of Si with 
respect to Fe. Results from simulations will be compared to the 
most recent observational data from the Chandra, XMM-Newton 
and Suzaku satellites. All the abunda nce values will be scaled to 
the the solar abundances provided by Grevesse & Sauvall (1998). 

To qualitatively appreciate the effect of different feedback 
mechanisms on the ICM enrichment pattern, we show in Figure 
|7] the maps of the emission-weighted Fe abundance for the dif- 
ferent runs of the g51 massive cluster. Here and in the following, 
we will rely on em ission-weighted estimates of metal abundances. 
Rasia et al. (2008) have shown that this emission-weighted estima- 
tor actually reproduces quite closely the values obtained by fitting 
the X-ray spectra of simulated clusters, for both Iron and Silicon. 
As for Oxygen, the emission-weighted estimator has been shown 
to seriously overestimate the corresponding abundance, especially 
for hot (T£, 3 keV) systems. 

In the run with no AGN feedback (upper panels of Fig|7}, we 
clearly note that including galactic winds produces a level of en- 
richment which is lower than that in the NF run outside the core re- 
gion, as a consequence of the lower level of star formation. There- 
fore, although galactic ejecta are known to be rather efficient in 
spre ading metals in the intergalac t ic medium at high re dshift {zlZ 2; 
e.g.. lOppenheimer & D ave 2008; ITescari et alj|2009l . Tornatore et 
al. 2009, in preparation), their effect is not strong enough to com- 
pensate the reduction of star formation within cluster regions. In the 
NF run we note the presence of highly enriched gas clumps which 
coincide with the halos of galaxies where intense star formation 
takes place. 

A rather different enrichment pattern is provided by AGN 
feedback (lower panels of Fig(7J. Despite the total amount of stars 
produced in these two runs is smaller than for the run including 
galactic winds, AGN feedback is highly efficient in spreading met- 



als at high redshift, mostly in correspondence of the peak of BH 
accretion activity. This demonstrates that AGN feedback provides 
a rather high level of diffuse enrichment in the outskirts of galaxy 
clusters and in the inter-galactic medium surrounding them at low 
redshift. 

4.1 Profiles of Iron abundance 

We show in Figure [8] the emission- weighted Iron abundance pro- 
files obtained by averaging over the four simulated clusters with 
T s i > 3 keV (left panel) and the five galaxy groups with T s i < 3 
keV (right panel), compared with observational results. Each panel 
reports the results for the four adopted feedback schemes. For rea- 
sons of clarity we report the la scatter computed over the ensemble 
of simulated clusters only for the AGN1 runs. 

As for rich clusters, sim ulation predictions are compa red with 
the observational results by iLeccardi & Molendil J2008al) . These 
authors analysed about 50 clusters with 3 keV, that were se- 
lected from the XMM-Newton archive in the redshift range 0.1 ^ 
z ^ 0.3. After carrying out a detailed modelling of the background 
emission, they recovered metallicity profiles out to ~ 0.4i?iso. 
The results of this analysis show a central peak of Zfb, followed 
by a decline out to 0.2i?i8o, while beyond that radius profiles are 
consistent with being flat, with Zpe — 0.3ZF e,p> using the solar 
abundance value by Grevesse & Sauval ( 1998) (~ 0.2 in units of 
the solar abundance by [Anders & Grevesse! jl989h , as reported by 
ILeccardi & Molendill2008ah . 

All our simulations predict the presence of abundance gradi- 
ents in the central regions, whose shape is in reasonable agreement 
with the observed one, at least for 0. 1-Riso ■ The lowest enrich- 
ment level is actually found for the NF run, despite the fact that this 
model produces the most massive BCGs. The reason for this lies 
in the highly efficient cooling that selectively removes the most en- 
riched gas, which has the shortest cooling time, thus leaving metal 
poorer gas in the diffuse phase. Runs with galactic winds (W) and 
AGN feedback (AGN1 and AGN2) are instead able to better reg- 
ulate gas cooling in central region, thus allowing more metal-rich 
gas to survive in the hot phase. For this reason, W and AGN runs 
predict profiles of Zf c which are steeper than for the NF runs in 
the central regions, 0.1i?iso- Quite interestingly, the effect that 
different feedback mechanisms have in displacing enriched gas and 
regulating star formation almost balance each other in the central 
cluster regions, thus producing similar profiles. However, the dif- 
ferent nature of SN-powered winds and AGN feedback leaves a 
clear imprint at larger radii. 

As for the runs with no feedback (NF), they produce a rather 
high level of enrichment out to ~ 0.3-Riso, while rapidly declin- 
ing at larger radii. In this model, the high level of star formation 
provides a strong enrichment of the gas in the halo of galaxies 
which will merge in the clusters. During merging, this gas is ram- 
pressure stripped, thus contributing to enhance the enrichment level 
of the ICM. The situation is different for the runs with winds. As 
already mentioned, galactic outflows are efficient in displacing gas 
from galactic halos at relatively high redshift, zjj 2, when they 
provide an important contrib ution to the enrichment of t he inter- 
galactic medium (IGM; e.g.. lOppenheimer & Dave 200j|). At the 
same time, this feedback is not efficient to quench cooling of en- 
riched gas at low redshift. As a consequence, no much enriched 
gas is left to be stripped by the hot cluster atmosphere from the ha- 
los of merging galaxies, thus explaining the lower enrichment level 
beyond O.l-Riso- 

As for the runs with AGN feedback, they produce a shape on 
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Figure 7. Maps of emission weighted Fe abundance in the g5 1 cluster for the runs without feedback (NF, top left), with winds (W, top right) and with AGNs 
(AGN1 and AGN 2, bottom left and b ottom right, respectively). Each map has a side of 2R v i r . Abundance values are expressed in units of the solar value, as 
reported by Grevesse & Sauval 1 1998), with color coding specified in the right bar. 




Figure 8. Comparison between the observed and the simulated profile s of emission-weighted I ron metallicity. Left panel: average Zp c profiles for galaxy 
clusters with T500 > 3 keV. Observational data points are taken from Leccardi & Molendi (2008a). Right panel: average Zp c profiles for the five simulated 
galaxy groups with T500 < 3 keV. Observational data points are taken from lRasmu^serT&Pomnar] j2007l) . In both panels different lines correspond to the 
average profiles computed for the different runs: no feedback (NF, blue short dashed), galactic winds (W, green long dashed), standard AGN feedback (AGN1, 
red dot-dashed), modified AGN feedback (AGN2, cyan solid). For reasons of clarity, we show with lcr errorbars over the ensemble of simulated clusters only 
for the AGN1 runs. 
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the abundance profiles quite close to the observed ones, with a flat- 
tening beyond ~ 0.27?iso- In this case, the effect of AGN feedback 
is that of displacing large amounts of enriched gas from star form- 
ing regions at high redshift (see also lBhattacharva et all2008r) and, 
at the same time, to efficiently suppress cooling at low redshift. The 
fact that the level of Zfc is almost constant out to Riso and beyond, 
witnesses that the main mechanism responsible for enrichment in 
this case is not ram-pressure stripping, whose efficiency should de- 
cline with cluster-centric radius. Instead, enrichment is this case is 
dominated by the diffuse accretion of pre-enriched IGM. We note 
that the AGN2 scheme tends to predict slightly higher Zp e values 
that AGN1. This is due to the effect of the more efficient radio- 
mode feedback, included in the former scheme, which provides a 
more efficient removal of gas from the halos of massive galaxies. 

Although models with AGN feedback produce the correct 
shape of the Iron abundance profiles, their normalisation is gen- 
erally higher than for the observed ones. This overproduction of 
Iron could be due to the uncertain knowledge of a number of in- 
gredients entering in the chemical evolution model implemented 
in the simulation code. For instance, differences between differ- 
ent sets of stellar yields turn into significant differences in the re- 
sultin g enrichment level (e.g jTornatore et a l. 2007; W iersma et al.l 
2009). Furthermore, a reduction of the Iron abundance can also be 
achieved by decreasing the fra ction of binary syst ems, which are 
the progenitors of SNe-Ia (e.g. jFabian et alj|2008f ). For these rea- 
sons, we believe that the shape of the abundance profiles, instead of 
their amplitude, should be considered as the relevant observational 
information to be used to study the impact that different feedback 
mechanisms have on the ICM enrichment pattern. 

In the right panel of Fig. [8] we compare the average 
Z-pc profiles of the five simulated galaxy groups with observational 
results fr om the analysis of 15 nearby g alaxy groups observed with 
Chandra iRasmu Tsen & Ponmanll2007l) . Also in this case, the pro- 
files from simulations have a slope quite similar to the observed 
one out to R ~ 0.3f?soo, although with a higher normalisation. At 
larger radii, the effect of AGN feedback is again that of providing 
rather flat profiles. This result is at variance with the observed pro- 
files in the outermost radii. Indeed, differently from rich clusters, 
galaxy groups apparently show a negative gradient of Iron abun- 
dance out to the largest radii covered by observations, with no ev- 
idence of flattening. This result further demonstrates the relevance 
of pushing observational determinations of the ICM enrichment out 
to the large radii, which it is mostly sensitive to the nature of the 
feedback mechanism. If confirmed by future observations, this re- 
sult may indicate that AGN feedback needs to be mitigated at the 
scale of galaxy groups, for it not to displace too large amounts of 
enriched gas. 

As pointed out bv lRebusco et al j |2005[). the central Z f,. peak 
in cool-core (CC) clusters (see also lDe Grandi et al.|[2004l) should 
be closely related to the star (light) distribution of the BCG. On 
the other hand, differences between the stellar mass profile and the 
metal abundance profile should be the signature of dynamic pro- 
cess es which mix and tra nsport metals outside the BCG. In this 
way, iRebusco et al] (2005) derived the amount of diffusion, to be 
ascribed to stochastic gas motions, which is required to explain 
the shall ower profiles o f Z_F e with respect to the BCG luminosity 
profiles. Roedig er et al.l J2007I) carried out simulations of isolated 
cluster-sized halos in which bubbles of high-entropy gas are in- 
jected to mimic the effect of AGN feedback. They showed that the 
gas diffusion associated to the buoyancy of such bubbles is indeed 
able to considerably soften an initially steep metallicity profile. 

In order to verify to what extent gas dynamical processes, as- 



sociated either to the hierarchical cluster build-up or to feedback 
energy release, are able to diffuse metals, we compare in Figure [9] 
the average Zp a profiles and the stellar mass profiles for the mas- 
sive clusters. As for the NF and W runs, we clearly see a central 
peak in the stellar mass density profile, related to the BCG, fol- 
lowed by a more gentle decline, which t races a diffuse halo o f intra- 
cluster stars surrounding the BCG (e.g., Murante et al .120071) . Quite 
clearly, the Zp c profiles are flatter than the distribution of stars at 
Since no feedback processes are at work in the NF 
runs, the flatter Zp a profile is due to the effect of selective cooling 
of highly enriched gas, rather than to stochastic gas motions. The 
same argument can be applied also to the runs including galactic 
winds (W), for which the kinetic energy provided galactic ejecta is 
not enough to fully regulate star formation in central cluster regions 
(see also Fig. [5}- At intermediate radii, 7? ~ (0.1 — 0.5)R V i r , pro- 
files of Iron abundance and stellar mass have quite similar slopes 
in both the NF and W runs. At even larger radii, instead, the 
Z-p a profile becomes again shallower than the stellar mass profile, 
a trend which persists even beyond the virial radius. At such large 
radii the shallower slope of Zp c can not be explained by the effect 
of gas cooling. It is rather due to the effect of gas dynamical pro- 
cesses which help mixing the enriched gas. 

As for AGN feedback, its effect is instead that of providing a 
similar slope for Z-p c and stellar mass profiles at R^z 0.2R V i r . This 
similarity is due to two concurrent effects. The first one is the sup- 
pression of star formation, which allows now central enriched gas 
to be pressurised by the AGN feedback, so as to leave a relatively 
larger amount of metal-enriched gas in the hot phase. The second 
one is the shape of the stellar mass profile, which is less concen- 
trated in the presence of AGN feedback. While the two profiles are 
similar in the central regions, the Zp a profile is instead much shal- 
lower at larger radii, _fix 0.3R V i r . Even an efficient AGN feedback 
can not be able to displace at low redshift significant amounts of 
enriched gas from the central regions of massive clusters beyond 
R V ir- Therefore, the flatter abundance profiles in the outskirts of 
rich clusters can only be justified by the action played by AGN 
feedback at high redshift. In fact, at Z^ 2 AGN efficiently removed 
enriched gas from galaxies and quenched star formation, thus pre- 
venting metals from being locked back in the stellar phase. 



4.2 When was the ICM enriched? 

In this section we discuss how the different feedback mechanisms 
change the cosmic epoch at which ICM was enriched. To this pur- 
pose, we define the average age of enrichment of a gas particle at 
redshift z as 



where the sum is taken over all time-steps performed by the simu- 
lation until redshift z, Amz,i is the mass in metals received by the 
particle at the i-th time-step, ti is the cosmic time of that time-step 
and mz(z) is the total metal mass received by the particle before 
z. According to this definition, a large value of enrichment age, 
at a given redshift, corresponds to more recent enrichment, while 
smaller values of i(z) indicate more pristine enrichment. In the 
limit in which all the metals are received by a particle at the consid- 
ered redshift z, then the enrichment age coincides with the cosmic 
age at z. Once computed for each gas particle, we then compute the 
mass-weighted mean of these ages of enrichment taken over all the 
particles having non-vanishing metallicity. 
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Figure 9. Comparison between the stellar mass density profiles, normalised to the corresponding value computed within -R200 (black short dashed lines), and 
the Iron abundance profiles (coloured lines). For reasons of clarity stellar profiles are shifted by a factor 10. All curves refer to the average profiles computed 
over the four massive clusters. From top left to bottom right panels we show results for the runs with no feedback (NF, blue short dashed), with galactic winds 
(W, green long dashed), with standard AGN feedback (AGN1, red solid) and modified AGN feedback (AGN2, cyan dot-dashed). 



Figure[T0]shows how the enrichment age of the ICM changes 
with the cluster-centric distance for the different feedback mod- 
els, for both rich clusters (top panel) and for poor clusters (bottom 
panel). In all cases, we note a decline towards small radii, although 
with different slopes. This demonstrates that gas in central regions 
has been generally enriched more recently than in the outskirts. In 
the innermost regions, the typical age of enrichment correspond to 
a look-back time of about 2-3 Gyr, while increasing to 5-8 Gyr 
around the virial radius. This confirms that metal enrichment in the 
central cluster regions receives a relatively larger contribution by 
star formation taking place in the BCG, and by stripping of en- 
riched gas from infalling galaxies, whose star formation has been 
"strangulated" only recently by the action of the hot cluster atmo- 
sphere. On the contrary, enrichment in the outskirts has a relatively 
larger contribution from high-z star formation, which provides a 
more widespread IGM enrichment. 

As for massive clusters, we note that the effect of galactic 
winds (W runs) is that of providing an earlier enrichment with re- 
spect to the model with no efficient feedback (NF runs), with a 
difference of about 1 Gyr, at all radii. This is quite expected, ow- 
ing to the efficient action that winds play in displacing metals from 
star forming regions at high redshift. The difference between W 
and NF runs is larger for the less massive clusters, consistent with 
the expectation that winds are more efficient in transporting metals 
outside shallower potential wells. 



Quite interestingly, AGN feedback provides a steeper radial 
dependence of the enrichment age, with rather similar results for 
the two alternative schemes of implementation (AGN1 and AGN2). 
In the central regions it is comparable to that of the NF runs, 
while becoming larger than that of the W runs for R> 0.3R vtr 
(^ O.GRvir) for rich (poor) clusters. The relatively recent enrich- 
ment age in central regions may look like a paradox, owing to the 
effect that AGN feedback has in suppressing low redshift star for- 
mation within the BCG (see Fig[T). However, we should remind 
that an aside effect of suppressing star formation is also that of 
preventing recently enriched gas from being locked back in stars. 
Therefore, although metal production in central regions is sup- 
pressed in the presence of AGN feedback, this effect is compen- 
sated by the suppression of cooling of recently enriched gas. The 
effect of AGN feedback in providing a more pristine metal en- 
richment becomes apparent in the outskirts of clusters. Indeed, the 
IGM in these regions has been efficiently enriched at high redshift, 
z ~ 3^1, when BH accretion reached its maximum activity level 
and displaced enriched gas from star forming galaxies. After this 
epoch, little enrichment took place as a consequence of the rapidly 
declining star formation, thus justifying the older enrichment age. 
Consistently with this picture, we also note that the increase of 
enrichment age with radius is more apparent for poorer clusters, 
where AGN feedback acted in a more efficient way. 
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Figure 11. Comparison between observed (solid squares with errorbars) and simulated (coloured symbols) relation between global Iron abundance and 
temperature. The left panel shows the results for the NF (dark blue) and W (light green) runs, while the right panel shows the results for the AGN1 (dark 
red) and AGN2 (light cyan) runs. In all cases, filled and open circles refer to the main halos of the resimulated Lagrangian regions with massive clusters and 
to the satellites, respectively. The main-halo galaxy groups are instead plotted with filled squares. Emission-weighted Iron abundance and spectroscopic-like 
temperature for simulatio ns are both computed within -Rsnn. Observational results refer to the sample of clusters observed with ASCA and analysed by 
Baumgartner et al. (2005). 
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Figure 10. Average age of enrichment as a function of the cluster-centric 
distance. The y-axis is for the look-back time at which the ICM was en- 
riched (in Gyr). The top panel is for the four massive clusters with M200 > 
10 15 /i -1 Mq, while the bottom panel is for the five galaxy groups with 
M200 — 10 14 h~ 1 M(7 ) . Different lines corresponds to the average age of 
enrichment computed within the set of simulated cluster: no feedback (NF, 
short dashed), galactic winds (W, long dashed), standard AGN feedback 
(AGN1, dot-dashed), modified AGN feedback (AGN2, solid). 



4.3 The metallicity - temperature relation 

In FigurefTTIwe present results on the abundance of Iron as a func- 
tion of cluster temperature, by comparing simulation results for 
the different feedback schemes with observation al results fro m the 
analysis of the ASCA Cluster Catalog ue (ACC: lHornell200lh car- 
ried out by Baumgart ner et al] ((2005). Due to the relatively poor 
angular resolution of the ASCA satellite, the extraction region for 
each cluster was selected to contain as much flux as possible. Since 
a unique extraction radius is not defined for the observed catalogue, 
we adopt -R500 as a common extraction radius to compute Iron 
abundances and temperatures of simulated clusters. We verified 
that adopting instead a larger extraction radius (e.g. R V i r ) slightly 
lowers the spectroscopic-like temperatures with out changing sub- 
stantia lly the results on Zy c ■ In their analysis, iBaumgartner et all 
(2005) partitioned their large sample of clusters in temperature 
bins. Then, a global value of Zp a was computed by combining all 
the clusters belonging to the same temperature bin. In this way, er- 
rorbars associated to the observational data points shown in Fig. 
QT]only account for the statistical uncertainties in the spectral fit- 
ting procedure after combining all clusters belonging to the same 
temperature bin, while they do not include any intrinsic scatter 
(i.e., cluster-by-cluster variation) in the Zf b -T relation. This has 
to be taken in mind when comparing them to simulation results, for 
which we did not make any binning in temperature. 

In both panels we see that simulations of hot (T > 6 keV) sys- 
tems produce values of Zf c which are above the enrichment level 
found in observations, ~ O.SZq. This result is consistent with 
the fact that Zp c profiles for simulated clusters have a larger nor- 
malisation than the observed ones (see Fig. [Hi. NF and W runs 
predict Zf & ~ (0.5 — 0.6)Z Fe q, with a slightly lower value, 
Z Fe ~ (0.4 - O.5)Z Fei for the runs with AGN feedback. 
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As for the simulations of the poorer systems, with T ~ (1 — 3) 
keV, they have instead a larger scatter. This indicates more diver- 
sity in the impact that gas-dynamical and feedback processes have 
in determining the enrichment pattern within smaller systems. The 
effect of AGN feedback for these systems is that of decreasing the 
value of Z-p c by about Q.IZfb.q- The rather large spread of abun- 
dance values and the limited number of simulated systems do not 
allow to establish whether simulations reproduce the increase of 
Zp e with ICM temperature for systems in the range T ~ (1 — 3) 
keV. 

Observations suggests that Zp c is almost independent of tem- 
perature above 6 keV, a trend which is in fact reproduced by simula- 
tion results. However, between 3 and 6 keV observed clusters show 
a drop in metallicity by about a factor two. Although simulations 
predict somewhat higher metallicity values at low temperatures, 
still it is not clear whether they reproduce the decrease, by more 
than a factor 2, suggested by observations. A potential complica- 
tions in comparing observational and simulation results for systems 
with T ~ 3 keV is that the spectroscopic value of Zf c for these sys- 
tems i s cont ributed by both K and L lines. As originally noted by 
iBuotej feOOOt) in the analysis of ASCA data, fitting with a single- 
temperature model a plasma characterised by a multi-temperature 
structure, with the colder component below 1 keV, leads to an un- 
derestimate of the Iron abundance (the so-calle d iron-bias; see also 
iMolendi & Gastaldellol200lMBuote et alj2003h . In this case, a sim- 
ple emission-weighted definition of Zp c from simulations may not 
be fully adequate. A correct procedure would require extracting a 
mock X-ray spectrum from simulated clusters, to be fitted with 
a multi-temperat ure (and multi-metallicity) plasma model (e.g., 
Rasia et al. 2008). Finally, one should also note that more recent 
determinations of the Z-Fe-T relation from XMM-Newton data, 
although based on a much smaller number of clusters, suggests 
a less pronounced decrease for systems hotter than 3 keV (e.g., 
IWerner et al1l2008h . There is no doubt that a systematic analysis 
of nearby clusters within the Chandra and XMM-Newton archives 
would help to confirm or disprove the metallicity-temperature re- 
lation based on ASCA observations. 

4.4 The Zsi /Zf c relative abundance 

The relative abundance of elements produced in different pro- 
portions by different SN types is directly related to the shape 
and possible evolution of the initial mass function. Furthermore, 
studying how relative abundances change with cluster-centric dis- 
tance provides insights on the different timing of enrichment 
and on how different metals, produced over different time-scales, 
ar e mixed by gas-dynamical proce s ses. ASCA data an alysed 
bylLoe wenstein & Mushotzky ( 1996). iFukazawa et al.l dl998h and 
Finog uenov et alj ( l2000h originally suggested that cluster outskirts 
are predominantly enri ched by SNe-II. A similar resu lt was found 
more recently also by iRasmussen & PonmarJ J2007), who anal- 
ysed XMM data for poor clusters with TJi 3 keV. On the other 
hand, Suzaku observat i ons of low temperature clusters and groups 
dSato et alj200ll2009l a. [2009l b) show instead a rather flat profile of 
Zsi/ZFe out to large radii, ~ 0.37?^,-, thus implying that SNe-Ia 
and SNe-II should contribute in similar proportions to the enrich- 
ment at different radii. 

In this section, we focus our attention on the Z^/Zfc ratio 
for low temperature clusters. We show in Figure [72l the emission- 
weighted map of the Zsi/Zp c ratio for the four different feedback 
schemes applied to the g676 cluster. For the run with no feedback 
(NF), we note that Zsi/Zp E is generally quite patchy. It reaches 



higher values in correspondence of high-density star forming re- 
gions, a feature which is also shared in different proportions by the 
other runs. Indeed, the products of SNe-II are released over a short 
time scale, since they are synthesised by massive stars. As a result, 
their distribution tends to trace preferentially the distribution of 
star-forming regions. On the other hand, SNe-Ia release metals over 
a longer time-scale. In fact, these stars have time to leave star form- 
ing regions, as a consequence of the same dynamica l effects which 
generate a population of inter-galactic stars (e.g., iMurante et al.l 
2007: IZibetti et alj |2005). thereby providing a more widespread 
enrichment pattern. Therefore, our simulations predict that diffuse 
intra-cluster stars provide a significant con tribution to the enric h- 
ment of the intra-cluster medium (see also lTornatore et alT 2007), 

A comparison of the maps for the NF and W runs shows that 
they have comparable levels of Z^/Zfc within the central regions. 
In these regions we expect that ram-pressure stripping is the dom- 
inant process in rem oving gas from merging galaxies (see also 
Kapf ereret all [2007) and, therefore, in efficiently mixing in the 
ICM the nucleosynthetic products of different stellar populations. 
On the contrary, in the cluster outskirts winds have been much 
more efficient in removing freshly produced metals from galaxies 
during the cluster assembly, therefore providing a relatively more 
widespread Si distribution, with respect to the NF case in which no 
galactic outflows are included. 

The runs with AGN feedback exhibit a behaviour in the clus- 
ter outskirts which is qualitatively similar to that of the W run, al- 
though with slightly higher values of Z^/Zfc ■ This indicates that 
AGN feedback has a higher efficiency in mixing SN-Ia and SN-II 
products. As for the central regions, the truncation of recent star 
formation by AGN feedback would lead to the naive expectation 
that a relatively higher value of Iron abundance with respect to Sil- 
icon should be found. The results shown in the bottom panels of 
FigHHlead in fact to the opposite conclusion, with a marked in- 
crease of Zsi/ZF C in the core regions. The reason for this lies again 
in the effect of selective removal of metal-enriched gas associated 
to cooling. Total metallicity of the gas around the BCG is domi- 
nated by SNe-II products. Therefore, gas more enriched by SNe-II 
has a relatively shorter cooling time. As a consequence, suppression 
of cooling in the core regions by AGN feedback tends to increase 
the amount of SN-II products in the ICM, thereby justifying the 
increase of Zsi/Z-F e with respect to the runs not including AGN. 

This qualitative picture is also confirmed by the profiles of 
Zsi/Zvc , that are plotted in Figure [T3] In each panel, simula- 
tion results show the average profile, computed over the 5 main 
relatively poor clusters, for each feedback model. In the left 
pa nel simulation results a r e com pared with observational results 
by IRasmussen & Ponmanl {2007) from the analysis of 15 nearby 
galaxy groups observed with Chandra. Also shown with the two 
horizontal lines are the values of Zsi/Z^c produced by SNe-Ia and 
SNe-II for a simple stellar population of initial solar metallicity, 
having a Salpeter IMF, using the same sets of yields adopted in our 
simulations. 

Observational data show a rather flat profile of Z^/Zfb with 
a value close to solar at small radii, followed by a sudden in- 
crease beyond ~ 0.2i?5oo- Taken at face value, this result would 
imply that at ~ 0.5i?5oo the enrichment is mostly contributed by 
SNe-II, while a mix of different stellar populations is required in 
the cluster centre. Although this result is qualitatively similar to 
that suggested by the visual inspection of the Zsi/Zp c maps of 
Fig ] 121 simulation resu l ts ar e quantitatively different from those by 
Rasmussen & Ponmanl d2007l) . Indeed, in no case simulations pre- 
dict an increase of ZsilZFe beyond 0.27?soo. While the maps sug- 
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Figure 12. Maps of the emission-weighted Zgi/Zp e distribution for the runs of the g676 cluster without feedback (NF, top left), with winds (W, top right), 
with AGNs ( AGN1 and AGN2, bottom le ft and right, respectively). The side of each map is 2R v i r . Abundance values are expressed in units of the solar value, 
as reported by Grevesse & Sauval 1 1998), with color coding specified in the right bar. 



gest that such an increase is also expected in simulations, it takes 
place only for R > R500 , with the largest values reaching at most 
the solar one. Within the radial range covered by Chandra obser- 
vations, both NF and W runs show rather flat profiles. Also in this 
case, we are more interested in the slope of the Zsi/Zp c profile, 
rather than in its normalisation. Indeed, too low values of ZsJZf^ 
in simulations by about 0.3 (in solar units) can be be either due to 
the choice of the IMF or to adopted stellar yields. 

As for the runs with AGN feedback the profiles of ZstlZFe 
confirm the expectation gained from the maps of Fig[T2]for a rel- 
ative increase of the Si abundance in the central regions. As for 
the behaviour at large radii, R^ 0.27?soo, also AGN feedback does 
not predict the pronounced i ncrease see in the Chandra data by 
Ras mussen & Ponmar] J2007t) . 

As already mentioned, this observational evidence for an en- 
hancement of SN-II pr oducts in the out s kirts of grou ps is not con- 
firmed by Suzaku data JSato et alj|2lMl2009kl2009l b) and XMM- 
Newton (Silvano Molendi, private communication). As shown in 
the left panel of Fig[T2] Suzaku observations of Abell 262, NGC 
507 and AWM 7 indicate a profile of Zsi/Zp c which is consistent 
with being flat out to 0.3R V i r (corresponding to about 0.5i?soo)- 
While it is not the purpose of this work to address the reason for 
the different results coming from different satellites, we want to 
stress the great relevance of tracing the pattern of ICM enrichment 
out to the largest possible radii. Indeed, this is the regime where 
gas-dynamical processes related to the cosmological build-up of 
clusters, past history of star formation and nature of feedback pro- 



cesses regulating star formation, all play a role in determining the 
distribution of metals. 

In summary, the analysis of chemical enrichment in our sim- 
ulations of galaxy clusters confirms that the resulting metallicity 
distribution in the ICM is given by the interplay between gas cool- 
ing, which tends to preferentially remove more enriched gas from 
the hot phase, feedback processes, which displace gas from star 
forming regions and regulate star formation, and gas dynamical 
processes associated to the hierarchical build-up of galaxy clusters. 
In particular, feedback implemented through the action of galactic 
winds powered by SN explosions or through energy extracted from 
gas accretion onto super-massive BHs, leave distinct features on the 
resulting pattern and timing of ICM enrichment. 



5 CONCLUSIONS 

We presented the analysis of an extended set of cosmological hy- 
drodynamical simulations of galaxy clusters aimed at studying the 
different effects that stellar and AGN feedback have on the ther- 
mal and chemo-dynamical properties of the intra-cluster medium 
(ICM ). Using a version of the Tree-SPH GADGET-2code dSpringell 
l2005h . w hich also include s a detailed description of chemical en- 
richment (Tornat ore et alj|2~007l) . we carried out simulations of 16 
clusters identified within 9 Lagrangian re gions extracted fr om a 
lower-resolution parent cosmological box (Dolag et al. 200^). All 
cluster simulations of this set have been run using different pre- 
scriptions for the feedback: without including any efficient feed- 
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Figure 13. Comparison between observed and the simulated profiles of Silicon abundance relative to Iron, Zgi/Zp c . In each panel, different lines correspond 
to the average profiles computed over clusters having T500 < 3 keV for the different sets of runs: no feedback (NF, blue short dashed), galactic winds 
(W, green long dashed), standard AGN feedback (AGN1, red dot-dashed), modified AGN feedback (AGN2, cyan solid). For reasons of clarity, we show 
with errorbars the r.m.s. scatt er over the ensemble of simulated clusters only f or the AGN1 runs. Obse r vation al points refer to the Chandra data analysed by 
iRasmussen & Ponmanl (2007) (left panel) and to the Suzaku data analysed by Sato et al. (2008, 2009a, 2009b) (right panel). The two horizontal lines show 
the relative abundance from SNe-Ia and SNe-II, computed for a simple stellar population (SSP) having solar initial metallicity and based on the same Salpeter 
IMF and set of yields as used in our simulations. 



back (NF runs), including only the effect of galactic winds pow- 
ered by supernova (SN) feedback (W runs), and including two dif- 
ferent prescriptions of AGN feedback (AGN1 and AGN2) based on 
modelling gas accretion on su per-massive black holes (BHs) hosted 
withi n resolved galaxy halos (Springel et al. 2005; Di Matt eo et al.l 
l2005h. The AGN 1 sch eme exactly reproduces the original model by 
JSpringel et ail2 005) for the choice of the parameters determining 
the feedback efficiency and the way in which energy is distributed. 
As for the AGN2 scheme, it assumes the presence of a radiatively 
efficient "r adio mode" phase w hen BH accretion is in a quiescent 
stage (e.g., Sijacki et al. 2007), also distributing energy to the gas 
particles surrounding BHs in a more uniform way. 

The main results of our analysis can be summarised as fol- 
lows. 

(a) AGN feedback significantly quenches star formation rate 
(SFR) associated to the brightest cluster galaxies (BCGs) at zj; 4. 
At z = the SFR in the AGN1 and AGN2 models is reduced by 
about a factor six. For a massive cluster with Af2oo — 10 15 /i _1 Mq 
we find SFR(z = 0) ~ 7 0M(^yr~ , thus not far from current ob- 
servational estimates (e.g. iRaffertv et al1l2006h . Although the two 
variants of AGN feedback produce similar results on the star for- 
mation rate, the AGN2 model is more efficient is reducing gas ac- 
cretion onto BHs. For this scheme, the resulting masses of the BHs 
sitting at the centre of the BCGs at z — are reduced by a fac- 
tor 3-5 with respect to the AGN1 scheme. Furthermore, increasing 
the radio-mode feedback efficiency from ej = 0.2 to 0.8 further 
reduces the mass of the central BH by about a factor 2.5, while 
leaving the level of low-z star formation rate almost unaffected. 

(b) AGN feedback brings the Lx~T relation in closer agreement 
with observational results at the scale of poor clusters and groups, 
thus confirming results from previous simulations based on AGN 



feedback JPuchwein etafl f2008). However, this is obtained at the 
expense of increasing the ICM entropy in central regions of groups 
above the level indicated by o bservational results (e.g. JSun et al.l 
120091 : ISanderson et alj|2~009bh . This entropy excess generated in 
central group regions corresponds in turn to a too low value of the 
gas fraction. 

(c) AGN feedback reduces by 30-50 per cent the fraction of 
baryons converted into stars, fstar, within f?soo- Simulation re- 
sults agree well with the observed value of fstar at the scale of 
groups. However, for rich clusters the fraction of stars within -R500 
from simulation (~ 30-40 per cent) is larger than the observed one 
(~ 10 per cent). 

(d) AGN feedback is quite efficient in pressurising gas in the cen- 
tral regions of galaxy groups, thereby generating temperature pro- 
files which are in reasonable agreement with the observed ones. 
Despite this success at the scale of groups, temperature profiles in 
the core regions of massive clusters are still too steep, even after 
including AGN feedback. 

(e) The presence of AGN feedback generates a rather uniform and 
widespread pattern of metal enrichment in the outskirts of clusters. 
This is the consequence of the improved efficiency, with respect 
to the runs without BH feedback, to extract at high redshift highly 
enriched gas from star forming regions, and, therefore, to enhance 
metal circulation in the inter-galactic medium. 

(f) Radial profiles of Fe abundance are predicted to be too steep 
at R?^ 0.17?i8o in runs including stellar feedback. Their shape is 
in much better agreement with the observed ones when includ- 
ing AGN feedback. The overall emission-weighted level of enrich- 
ment within massive clusters is Zfb — 0.5 and 0.6 for runs with 
and without AGN feedback, respectively. Such values are generally 
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larger than those, Zfz — 0.3, reported from ASCA observations 
(Baumgartner et al. 20051). 

(g) The distribution of elements mostly produced by SNe-II over 
a relatively short time-scale is more clumpy than the Iron distri- 
bution, which has a larger contribution from SNe-Ia, that release 
metals over a long time-scale. This is interpreted as due to the ef- 
fect of enrichment from stars belonging to a diffuse intra-cluster 
population. Therefore, simulations predict that a sizable fraction of 
the ICM enrichment is p roduced by intra-cluster stars, in line with 
observational evidences dSivanandam et alj 2009). 

(h) The runs with no feedback (NF) and with galactic winds (W) 
predict similar patterns of Zsi/Zp c within Rsoo- Silicon abundance 
is enhanced in the outer regions by the action of galactic winds, 
thanks to their efficiency in transporting gas enriched by a elements 
from star forming regions. AGN feedback has the effect of increas- 
ing Zs[/Zf c as a consequence of the suppression of star formation, 
which would lock back in stars gas surrounding star forming re- 
gions, and of the efficient removal of enriched gas from galactic 
halos. 

(i) In no case we find that profiles of Zsi/Zp c have a rising 
trend beyond ~ 0.27?soo. No strong conclusion can be drawn 
from a comparison with data, owing to discrepant indications 
from different observational results on the radial dependence of 
Zs[/Zf c . Suppression of star formation with AGN feedback 
causes Zsi/Z-pe to increase at small radii, ^ O.I-R500, a feature 
which is not seen in observational data. This suggests that the im- 
plementation of AGN feedback in our simulations may not provide 
enough gas mixing the central regions of clusters and groups. 

Our analysis lend further support to the idea that a feed- 
back source associated to gas accretion onto super-massive BHs 
is required by the observ ational properties of the ICM (e.g. 
iMcNamara & Nu lsen 2007). However, our results also show that 
a number of discrepancies between observations and predictions of 
simulations still exist, especially within the core regions of mas- 
sive clusters. This requires that a more efficient way of extracting 
and/or thermalising energy released by AGN should be introduced 
in richer systems. A number of observational evidences exists that 
AGN should represent the engine which regulates the structure of 
core regions of clusters and groups. However, observations also 
provide circumstantial evidences that a number of complex physi- 
cal processes, such as injection of relativistic particles and of tur- 
bulence associated to jets, buoyancy of bubbles stabilised by mag- 
netic fields, viscous dissipation of their mechanical energy, thermal 
conduction, should all cooperate to make AGN feedback a self- 
regulated process. In view of this complexity, we consider it as 
quite encouraging that the relatively simple prescriptions for en- 
ergy thermalisation adopted in our simulations provide a significant 
improvement in reconciling numerical and observational results on 
the ICM thermo- and chemo-dynamical properties. 

Clearly, increasing numerical resolution thanks to the ever in- 
creasing supercomputing power would require including a proper 
description of the above processes. For instance, the current imple- 
mentation of AGN feedback neglects the effect of kinetic energy 
associated to jets. The typical scales of ~ 20 kpc at which kinetic 
feedback is expected to dominate within clusters (e.g. Popej |2009l) 
are only marginally resolved by hydrodynamics in our simulations, 
thus making the assumption of a purely thermal feedback a reason- 
able one. However, the improved numerical resolution expected to 
be reached in simulations of the next generation needs to be ac- 
companied by a suitable description of injection of jets, for them 



to provide a physically meaningful description of the interplay be- 
tween BH accretion and ICM properties. 

The results presented in this paper further demonstrate that 
different astrophysical feedback sources leave distinct signatures 
on the pattern of chemical enrichment of the ICM. These differ- 
ences are much more evident in the outskirts of galaxy clusters, 
which retain memory of the past efficiency that energy feedback 
had in displacing enriched gas from star-forming regions and in 
regulating star formation itself. However, characterisation of ther- 
mal and chemical properties in cluster external regions requires X- 
ray telescopes with large collecting area and an excellent control of 
the background. While Chandra, XMM and Suzaku will be pushed 
to their limits in these studies in the next few years, there is no 
doubt that a detailed knowledge of the ICM out the cluster virial 
boundaries has to awa it for the advent of the next generation of 
X-ray telescopes (e.g.. lGiacconi et al . 2009 : 1 Arnaud etal .1120091) . 
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